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Abstract. Nonlinear perturbation of Fuchsian systems are studied in a region 
including two singularities. It is proved that such systems are generally not analyt- 
ically equivalent to their linear part (they are not linearizable) and the obstructions 
are found constructively, as a countable set of numbers. Furthermore, assuming a 
polynomial character of the nonlinear part, it is shown that there exists a unique 
formal " correction" of the nonlinear part so that the " corrected" system is formally 
linearizable. 

Normal forms of these systems are found, providing also their classification. 



1. Introduction 

1.1. Setting. Nonlinear perturbations of Fuclisian systems are, in the present con- 
text, differential systems of the form 

(1) ^ = M(x)u + g(x,u) 

studied for u G C"' small, and x in a simply connected domain D C C which includes 
singular points of the matrix M{x). 

It is assumed that the linear part of ([T]) 

(2) ^ = M(x)w weC^, xeC 

dx 

is a Fuchsian system. Then that all its singularities in the complex plane, including 
the point x = oo, are of Fuchsian type (see the Appendix, §5.11 for details). 

The function g(x, u), which gathers the nonlinear terms, is assumed to be holo- 
morphic, or, it can have (at most) simple poles at the singularities of M[x): 

(3) g(x,u) = ^^f(x,u) 

where: Q{x) = ]^(x — Xj), with Xi (z D singular for M, f (x, .) has a zero of order 
two at u = and f is holomorphic for x E D and small u. 

The case when the domain D contains just one singular point of M (x) was studied 
in [1] and the main results are summarized in §1.3.21 
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The present paper considers the case when M{x) has three singularities in the 
extended complex plane. Placing one at oo, the other two are conventionally placed 
at X = 1 and x = —1 (but their location can be arbitrarily placed using a linear 
fractional change of the variable x). In this case the matrix M has the form 

(4) Mix) = ^A + A,BeM,iC) 

X — 1 X + 1 

and 

(5) Q{x) =x^ -I 

(the condition that a; = oo is also of Fuchsian type implies that the matrices A and 
B do not depend on x). 

The systems studied have therefore the form 

du f 1 , 1\ 

where f is analytic for x in the simply connected domain D 3 {±1} and for u G C^, 
|u| < r, f(x, ■) has a zero of order two at u = 0. 

In addition it is also assumed that the nonlinear part is of polynomial type in x, 
in the sense that 

(7) f(x, u) = fin(a;)u'" with fni(x) polynomials 

m|>2 

1.2. Motivation. The question of linearization and, more generally, of equivalence, 
is a fundamental problem in the theory of differential equations. 

Besides the clear theoretical interest, linearization and equivalence are used as 
instrumental methods in control theory, and in devising algorithms for numerical 
and symbolic calculations (see, to cite just a few authors, [2]- [6]). Often the meth- 
ods used are developed from the method of equivalence introduced by Cartan [7] to 
decide whether two differential structures can be mapped one to another by a trans- 
formation taken in a given pseudogroup [8j, In the case of differential equations, 
the method was used for regular systems and extended for a neighborhood of one 
singular manifold [TU], [TT] . 

The present paper uses direct and constructive methods in the study of lineariz- 
ability and normal forms for equations in the class ([6]) in domains which include two 
singular points. 

The study of vector fields with an eigenvalue equal to 1 at a singular point can be 
reduced to the study of a Fuchsian system near one singularity ([8]) by ehminating 
time in a time-independent system. Similarly, the study of vector fields in regions 
containing two, or more singular points is reducible to that of equations of type ([T]). 
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More generally, the study of Hamiltonian systems with polynomial potentials near 
particular, periodic, or doubly-periodic solutions can be reduced to the study of 
nonlinear perturbations of Fuchsian systems [12] , [13] . 

The result of Theorem[3] has a very interesting similarity with needed corrections 
found in other problems: Ecalle and Vallet showed that resonant systems are lin- 
earizable after appropriate correction [14J; also Gallavotti showed that there exists 
appropriate corrections of Hamiltonian systems so that the new system is integrable 
[T5] . convergence being proved later by Eliasson [TB]. This suggests the possible 
existence of an underlying general structure. 

1.3. Linearization in a neighborhood of one regular singular point and 
connection to integrability. 

1.3.1. Notations. For w = (wi, . . . , Wd) G C'^ and m G N*^ denote 

w™ = . . . w^'' , |m| = mi + . . . + rrid , m • w = miWi + . . . + irtdWd 

1.3.2. Region containing one regular singularity. Consider a nonlinear perturbation 
of a system with one regular singular point in the complex plane: 

(8) ^ = -L(x)u + if(x,u), ugC"^, xGC 
ax X X 

where L{x) is a matrix analytic at a; = 0. 

In nonresonant cases, after an analytic change of variables, it can be assumed that 

the matrix L{x) is constant (see §5.21 for details): 

(9) ^ = -Lu + -i(x,u), ueC'^,xeC 
ax X X 

f(x, u) is assumed analytic for small u: |u| < p and for x in a domain Dr which is 
either a disk centered at the origin: |x| < r', or an annulu£| r" < \x\ < r'. 

Such systems are analytically linearizable if a Diophantine condition is satisfied 

[IJ: 

Theorem 1. Assume that the eigenvalues fii, ...,fid of the matrix L satisfy the Dio- 
phantine condition: there exist C,^ > so that 



(10) 



n ■ /X + / — yU< 



> C(|n| + 



for all I & N if Dr is a disk (respectively all I & Z if D^ is an annulus), and for all 
s G {1, d], and n G N'^ with |n| > 2. 

Then the system ^ is analytically equivalent to its linear part w' = for x E Df 
(a sub-disk, respectively, sub-annulus, of Dr: < f < r) and |u| < p < p. 



^Allowing the nonlinear part to be singular at a; = accommodates systems corresponding to 
higher order equations. 
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Remark 1. 

(i) The condition ^^) is satisfied by generic matrices L in the sense that the 
Lebesgue measure of the set of points G which do not satisfy a condition [T^) 
is zero, see ^7\,Ch.6,^ 3. 

(a) The domain Df can be made arbitrarily close to Dj. for p small enough. 

(Hi) The analytic equivalence map with identity linear part is unique if no eigen- 
value fij is integer. 

1.3.3. Region containing two singularities of Fuchsian type. Consider equation ([6]). 
By Theorem [T], there (generically) exists an equivalence map of this equation to its 
hnear part ([2]), for short, a linearization map, analytic at u = and for x in a 
neighborhood of 1. This map is unique under the requirement of analyticity at x = 1 
(see Remark[T]). Similarly, there exists a unique linearization map which is analytic 
at X = —1. Generally these two maps do not coincide: the map which is analytic at 
X = ±1 is (usually) ramified at x = =f1. Indeed, the first term of the right-hand-side 
of fl24|) below displays the ramification. 

Therefore equations ([1]) are not necessarily linearizable in a region containing more 
than one singularity. 

Absence of linearizability has important consequences for the system, since it im- 
plies non-integrability (at least generically) [TH]: 

Theorem 2. Consider the nonlinear equation in one dimension: 



for X in a domain D containing x = ±1 and small u. 

If equation 07]) is not analytically linearizable then for generic a, h (precise condi- 
tions are given in [18jj solutions have dense branching and therefore no single-valued 
integrals exist. Among integrable cases, first integrals are not meromorphic (generi- 
cally). 

It is therefore important to find what are the linearization criteria for equations 
([1]), and furthermore to describe the equivalence classes - to find normal forms for 
these equations. 



2.1. Assumptions. Consider the system ([6]). The following conditions (a)-(c), sat- 
isfied by generic matrices, are assumed. 

(a) The eigenvalues of A, respectively B, satisfy the Diophantine condition fllOp . 
and (b) none of these eigenvalues is an integer. 





2. Main Results 
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(c) The eigenvalues Ai, . . . , of the matrix A + B satisfy the following nonreso- 
nance condition: 

(12) k + n-X-\j^O for all n G N'^ , A; G N , j = l,...,d 

2.2. Obstructions to linearization and existence of corrections. 

Theorem 3. Consider the system ^ under the assumptions of ^2.1[ 
Then there exists a unique "correction" <p{u) ofi{x,u.) as a formal series 

(13) 0(u) = 0^u- 

|m|>2 

SO that the "corrected" equation 

, N dW ( \ ^ 1 \ 1 r„/ N , / M 

is linearizable by a formal series 

(15) u = H(x,w) = w+ ^ hm(a;)w™ 

|m|>2 

where h^{x) are functions analytic on D. 
Moreover, hni(a;) are polynomials. 

Remarks 

1. The coefficients (f)^ G of the correction <j) (see ( IT3l) ) represent the obstruc- 
tions to linearization, in the sense that the system (jSj) is linearizable by a formal 
power series if and only if all these vectors are zero. 

Obviously, if no formal power series f|T5|) exists, then analytic equivalence is pre- 
cluded too. 

2. are obtained constructively, see fl20l) . (|2T|) . Lemmas[5] andO 

3. If the nonresonance condition f|T2|) does not hold then existence, or uniqueness, 
of corrections may fail, as illustrated in the Appendix §5.31 

The condition that f(-,u) be of polynomial type is expected to generalizable to 
analytic functions. In a forthcoming paper [L9l it is shown that if the matrices A and 
B are simultaneously diagonalizable, then Theorem^ holds more generally, for f (x, u) 
analytic in x. Furthermore, it is also shown in [12] that if in addition, the eigenvalues 
of the matrices A and B have positive real part, then all the series converge and the 
corrected system is analytically linearizable. 

2.3. Formal normal forms. Since equations ([6]) are not necessarily linearizable, 
then two such equations are not necessarily equivalent. The following theorem finds 
normal forms for equations 
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Theorem 4. Consider the system ^ under the assumptions of ^2.1[ 

Then there exists a unique formal series ip{w) = X]|m|>2 'V'm'^™ ■^^ '^^^^ ® 
equivalent to 

(16) w'=(^A+^b)w + ^V'(w) 

through a formal series 

(17) u = H(x,w) = w+ ^ h„(2;)w™ 

|m|>2 

with hm(a;) analytic functions on D. 
Moreover, hni(x) are polynomials. 



3. Proofs 

3.1. Proof of Theorem [3l 

3.1.1. The recursive system. A map of the form (|T5l) . u = H(2;, w) = w + h(x,w), 
is a linearization map of (EI) if and only if h satisfies the nonlinear partial differential 
equation 

(18) a^h + dwh Mw = Mh + — — [f (x, w + h) - 0(w + h)] 

where M is given by (jlj). 

Denote by h„ the homogeneous part of degree n of a function h(2;,w): 

(19) h„(x,w)= hm(x)w- 

|m|=n 

(and a similar notation is used for f (x, w) etc.). 

Writing solutions of (fT8|) as power series in w in we obtain h„ recursively in n. 
The functions {hm}|jn|=n satisfy the linear differential system 

(20) d^K + d^K Mw - Mhn = R„(x, w) 

x^ — 1 

where 

(21) R„ = f„ - </)„ + R„ 
with R„ a polynomial in (p^, h^, f^ with |m| < n. 
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3.1.2. Numerical obstructions to linearizability: an analytic illustration. The exis- 
tence of the correction 0(u) is estabhshed below in §3.1.31 using algebraic methods 
which rely on the polynomiality assumption ([7]). It is interesting, however, to look at 
the analytic context, as a hint on why existence, and uniqueness, of this correction 
is to be expected. 

Denote by Y{x) a fundamental matrix of the linear system ([2]): Y' = MY. If 
Q^^Y^^ is Lebesgue integrable at x = ±1, then a particular solution of (1201) is given 
by 

(22) K{x,w) = Y{x)J Q{ty^Y{t)-^R^{t,Y{t)Y{xy'w)dt 

Let G be the monodromy matrix of ([2]) at x = —1: after analytic continuation 
along a closed loop around x = —1 the matrix Y{x) becomes y4C(_i)F(x) = Y{x)G. 
The analytic continuation of (l22l) on a closed loop around x = —1 yields 

(23) Aq_i)h„(x,w) = r(x)G J Q{t)-^G-^Y{t)-^Rn{t,Y{t)GG-^Y{x)-^w)dt 

= h„(a;, w) 

which means that fl22!) are the homogeneous terms of the unique linearization map 
of ([6]) which is analytic at x = — 1. 
Rewriting fl22l) as 

(24) K{x,w) = Y{x)j Q-^Y-^Rndt + Y{x) Q-^Y-^R^dt 



a similar argument shows that the last term in fl24l) is the homogeneous term of 
the unique linearization map which is analytic at x = 1. Then u = w + h(x, w) is 
analytic at both x = 1 and x = — 1 if and only if 

(25) j Q{t)~^Y{t)-^Rn{t, Y{t)Y{x)-^w) dt = for all w G n>2 

The recursive formulas ( 125|) contain the obstructions to linearization: there are d 
numerical conditions (one vector in C'') for every m G N'^, |m| > 2. 



3.1.3. Existence of the correction 0. The existence of (as a formal series) is es- 
tablished directly on the differential equation (I2U]) with (12 ip . rather than solving 
fl25|) . This result, which completes the proof of Theorem[3l follows directly from the 
following lemma: 

Lemma 5. Denote by Vn the space of vector-valued polynomials in w, homogeneous 
degree n (Vn C C'^[wi, . . . , Wd]). 
Assume /TUL) holds. 
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Then for any F = F(a;, w) G Vn[x\, there exists a unique (f) E Vn so that the 
differential equation 

(26) d^P{x, w) + dwP(x, w) M{x)w - M{x)P{x, w) = ^ [F(x, w) - 0(w)] 

has a solution P G 

The proof of Lemma E] is given in §3.1. 4[ 

Note that if F has degree + 1 in x, then the solution P of fl2Bl) has degree k. 



Remark 2. The only solutions of l\2b]) which are analytic at both x = 1 and x = —1 
are polynomials, and they are given by the following generalized Rodrigues formulas. 
Consider the polynomials Fk+i given by 

jk+i 

Ffc+i (x,w;q) = Q{x)Y{x)^ [Q{t)''Y{t)-'q{Y{t)Y{x)-'w)] , q g P„ , G N 
(these, together with Vn, generate the vector space Vn[x\). Then 

(27) P.(^,w;q) = y{x)^[^^ [Q{tfY{tr\{Y{t)Y{x)-'^)\ 

is the polynomial solution of / f^) for F = F^+i and <p = 0. 

The functions [2l\ ) are indeed polynomials in x, of degree k. For example 

Po (x, w; q) = q(w), Pi (x, w; q) = {Q' - QM) q(w) + Q rfq(w) Mw 
P2 (x, w; q) = [{Qy + Q\M^ - M') - 2{Q^y M] q(w)+2 [(Q^)' + q'^m] rfq(w) Mw 

+Q2 c/q(w) (M' + M2)w + g2 dP^^iy^) (Mw, Mw) 
3.1.4. Proof of LemmalM The proof uses the following solvabihty Lemma: 

Lemma 6. Let A be a d x d matrix with eigenvalues satisfying the nonresonance 
condition fT^) . Let J\ be the following linear operator on Vn.: 

(28) ( Jap) (w) = (dp) Aw - Ap 
Let I G N. 

(i) If A is diagonal, then the operator I + is diagonal, with eigenvectors w^ej, 
and corresponding eigenvalues I + m ■ \ — Xj , for all j = 1, . . . ,d and m, |m| = n. 

(a) If A is in Jordan normal form, then I + Ja is in Jordan normal form and the 
eigenvalues are as at (i). 

(a) As a consequence, the operator I + Ja is one-to-one and onto Vn- 

Lemma [6] follows by calculation (see also |20] Ch.5, §22). 

To prove Lemma \5\ denote by L and N the following linear operators on Vn' 

L = J(A+B) and = J(a-b) (see 
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Let F = F(x, w) G Vn[x] be a polynomial in x degree k + hF = with 
fj G Vn- Identifying the coefficients of x, a solution P of (1261) has degree k and, if 
P = X]j=o^"'Pi then must satisfy the recursive system 

{k + L)pk = ffc+i 
/^QN {k-l + L)pk-i = Npk + ffe 

^ ^ (/ - 1 + L)pi^, = Npi + + l)pi+i + fi for 1 < / < A; - 1 

= A^po + pi + fo - 

By Lemma [6] the system (!29|) can be solved recursively for p^, then p^-i, and all 
the way to po, for any ffc+i, ffc, . . ., fi. The last relation in (l29l) determines uniquely 

0. □ 

3.2. Proof of Theorem m The map (fT7|) is an equivalence map of (P) and ([TB|1 iff 
h satisfies the nonlinear PDE 

(30) d^h + d^h Mw - Mh = — [f (x, w + h) - V(w) - (d^h) V'] 

x'^ — 1 

where M is given by (jlj). As in the proof of Theorem [3l expanding in power series 
in w and denoting by h„ the homogeneous part of degree n in ( JT5l) we obtain h„ 
recursively as solutions of the equations 

(31) d,K + d^hnMw - MK = ^^Rn{x,w) {n>2) 
where 

(32) R„ = f„ - + Rn 

with R„ a polynomial in '0^, hm, fm with |m| < n. 

The terms i/)^ are recursively determined from (I3T!) . (l32l) using Lemma [5] for F = 

fn + Rn- n 
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5. Appendices 

5.L Singular points of Fuchsian type of a linear differential equation. Con- 
sider a linear differential equation 

(33) ^ = M(x)w w G xeC 

dx 

A point a; = Xo is called a regular point of ( l33l) if the matrix M{x) is analytic at 

X = Xq. 
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A point X = Xq is called singular of Fuclisian type if M{x) has an isolated singu- 
larity at X — Xq and there exists a constant k so that all solutions of the system, 
in every sector of the complex x-plane with vertex at Xq grow no faster than x'' as 
x^xo (see [TTj, Ch.7, §2.1). 

The equation (1551) is called Fuchsian if 

^ 1 
M{x) = J2 -^—^j 

j=l 3 

in which case solutions of fl33l) are analytic at all points except and oo and 

these points are of Fuchsian type. It turns out that in this case all solutions have 
formal asymptotic expansions as series in powers and logs, as x —>■ aj, or x —>■ oo 
|21] . Moreover, these series converge. 

For generic nonlinear perturbations of systems fl33l) near a singular point of Fuch- 
sian type Xq, solutions also have convergent asymptotic expansions in powers of x as 
X —>■ Xq hj Theorem[TJ For simplicity, the point x = Xq will be also called a singular 
point of Fuchsian type (of a nonlinear system). 

5.2. The matrix L{x) can be assumed constant in ([§]). Indeed, if the eigen- 
values of L{0) are nonresonant (in the sense that no two eigenvalues differ by an 
integer) then the linear part of ([8]): 

(34) 5 = -^(^)w 



diX X 



is analytically equivalent to 



(35) ^ = 1l(0)w 

for small x and w [T7|. Denoting w = H{x)w, the analytic equivalence map between 
and (155]) . the same map transforms ([S]) into a system with the same regularity, 
but with L constant: the map u = H{x)u is an analytic equivalence of ([H]) and 

(36) ^ = -L{0)u+-H-^~f{x,Hu) 

ax X X 

5.3. Examples in which condition (I12p does not hold. To illustrate that exis- 
tence of correction, or their uniqueness, may fail if the condition (fT2l) is not satisfied 
consider the simplest nonlinear equations: Riccati equations (in dimension one) 



(37) u'{x) = ( ^ + ^ ) n{x) + 4^ u\x) 

X— 1 X'^ 
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For n = 2 the equation (120|) is 

(38) /.'.(X) + + = 

\x — 1 X + 1 J X'^ — 1 

Take for simplicity f{x) = mx'^ + nx and a = b and note that for a = b = —1/2 
the condition f|T2|) is not satisfied. A simple calculation yields that equation fl38l) has 
an analytic solution for 02 = fn/ (2a + 1) and this solution is h2{x) = m/ (2a + l)x + 
n/{2a). Furthermore 0„ = for n> 2 and it can be easily checked that /i„ = K^^^ ■ 

These formulas certainly cannot be used in the resonant case a = 6 = — 1/2. In 
this case the general solution of (!38l) is 

(39) h2{x) = m [(1 — x^Y^"^ arcsinx — x] — n + (f)2X + C(l — x^)^^^ 

If m 7^ then no solution fl39l) is single- valued at both x = ±1, no matter what 
02 is: existence of corrections fails. 

By contrast, for m = uniqueness fails, since the solution ( l39l) with C = is 
analytic for any 02- 
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